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<N ■ Abstract 



This note describes a Mathematica interface for Fortran code generated by Form- 
er^ \ Calc. The interfacing code is set up automatically so that only minuscule changes in 
the driver files are required. The interface makes a function to compute the cross- 
section or decay rate available in Mathematica. This function depends on the model 
parameters chosen for interfacing in the Fortran code. 



O 

^! 1 Introduction 

^| The FeynArts/FormCalc system [H|2] is a fairly automated system for generating Fortran 
code to numerically evaluate a wide class of scattering and decay processes. Thus far, this 
Fortran code was a stand-alone, command-line driven application. Many users added extra 
functionality, either to the existing code, or by providing their own frontend code [3]. 

The new Mathematica interface described herein turns the stand-alone Fortran code 
into a Mathematica function for evaluating the cross-section or decay rate as a function of 
user-selected model parameters. The benefits of such a function are obvious, as the whole 
instrumentarium of Mathematica commands can be applied to them. For example, it is 
quite straightforward, using Mathematica's FindMinimum, to determine the minimum (or 
maximum) of the cross-section over a piece of parameter space. 

Interfacing is done using the MathLink protocol. The changes necessary to produce a 
MathLink executable are quite superficial and affect only the file run.F, where the user 
has to choose which model parameters are interfaced from Mathematica. 

To make the obvious even clearer, the cross-section is not evaluated in Mathematica, 
but in Fortran, and only the numerical results computed by the Fortran code are transferred 
back to Mathematica. One thing one cannot do thus is to increase the numerical precision 
of the calculation using Mathematica commands like SetPrecision. 

The Fortran results are delivered to and accessed in Mathematica in a way very similar 
to the ReadData utility which has been part of the FormCalc distribution for quite a while 
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already. The interfacing code for Mathematica is included from FormCalc Version 5.1 
onward. 

2 Setting up the Interface 

The model parameters are specified in the file run.F. Typical definitions for stand-alone 
code look like (here from an MSSM calculation): 

#define L00P1 do 1 TB = 5, 50, 5 
#define L00P2 MAO = 500 

These lines declare TB (= tan/3) to be scanned from 5 to 50 in steps of 5 and set MAO 
(= M A o) to 500 GeV. To be able to specify TB in Mathematica instead, the only change is 

#define L00P1 call MmaGetReal(TB) 

Such invocations of MmaGetReal and its companion subroutine MmaGetComplex serve two 
purposes. At compile time they determine with which arguments the Mathematica function 
is generated (for details see below), and at run time they actually transfer the function's 
arguments to the specified Fortran variables. 

• call MmaGetReal (r) 

read the real parameter r from Mathematica, 

• call MmaGetComplex (c) 

read the complex parameter c from Mathematica. 

Note that without a separate MmaGetReal call, MAO would still be fixed by the Fortran 
statement, i.e. not be accessible from Mathematica. 

Once the makefile detects the presence of these subroutines, it automatically generates 
interfacing code and compiles a MathLink executable. For a file run . F the corresponding 
MathLink executable is also called run, as in the stand-alone case. This file is not started 
from the command-line, but used in Mathematica as 

Install ["run"] 

3 The Interface Function in Mathematica 

After loading the MathLink executable with Install, a Mathematica function of the same 
name is available. For definiteness, we will call this function 'run' in the following since 
'run.F' is the default parameter file. This function has the arguments (grey parts are 
optional) 
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run[sqrtS, argl, arg2, options] 



run[{sqrtSf rom, sqrtSto, sqrtSstep}, argl, arg2, options] 

The first form computes a differential cross-section at y/s = sqrtS. The second form 
computes a total cross-section for energies \fs varying from sqrtSf rom to sqrtSto in steps 
of sqrtSstep. This is in one-to-one correspondence with the command-line invocation of 
the stand-alone executable. 

The argl, arg2, . . . , are the model parameters declared automatically by the presence 
of their MmaGet {Real, Complex} calls (see above). They appear in the argument list in the 
same order as the corresponding MmaGet{Real, Complex} calls. 

Other parameters are specified through the options (default values are written on the 
right-hand sides): 

• Polarizations -> "UUUU", the polarizations of the external particles, specified as 
in the stand-alone version, i.e. a string of characters for each external leg: 



• Serial -> O, a range of serial numbers, specified as {serialfrom, serialto, 
serialstep} (grey parts optional). The concept of serial numbers, used to parallelize 
parameter scans, is described in Ref. ^j. This option applies only to parameters 
scanned by Fortran do-loops in the parameter statements. Parameters read from 
Mathematica are unaffected by this option. 

• SetNumber -> 1, a set number beginning with which parameters and data are stored 
(see next Section). 

• ParaHead -> Para, the head under which parameters are stored, i.e. parameters are 
retrievable from parahead [setnumber] (see next Section). 

• DataHead -> Data, the head for the data storage, i.e. data are retrievable from 
datahead [setnumber] (see next Section). 

• LogFile -> "", a log-file to save screen output in. An empty string indicates no 
output redirection, i.e. the output will appear on screen. 

4 Return values, Storage of Data 

The return value of the generated function is an integer which records how many parameter 
and data sets were transferred. Assigning parameter and data sets as the data become 
available has several advantages: 



U unpolarized, 

T transversely polarized, 



L longitudinal polarization, 
- left-handed polarization, 
+ right-handed polarization. 
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• the return value of run is an integer rather than a large, bulky list, 

• the parameters corresponding to a particular data set are easy to identify, e.g. 
Para [4711] contains the parameters corresponding to Data [4711] , 

• most importantly, if the calculation is prematurely aborted, the parameters and data 
transferred so far are still accessible. 

Both, the starting set number and the heads of the parameter and data assignments can 
be chosen with the options SetNumber, ParaHead, and DataHead, respectively. 

The parameters which are actually returned are chosen by the user in the PRINTn 
statements in run.F in much the same way as parameters are selected for printout in the 
stand-alone code. To specify that TB and MAO be returned, one needs the definitions 

#define PRINT1 call MmaPutRealC'TB" , TB) 
#define PRINT2 call MmaPut Real ("MAO" , MAO) 

Notwithstanding, parameters can still be printed out, in which case they end up in the log 
file (or on screen, if no log file is chosen). To transfer e.g. TB to Mathematica and print it 
out, one would use 

#define PRINT1 call MmaPutRealC'TB", TB) 
#define PRINT2 SHOW "TB", TB 

An analogous subroutine exists of course for complex parameters, too: 

• call MmaPutReal(s, r) 

transfer the real parameter r to Mathematica under the name s, 

• call MmaPutComplex(s , c) 

transfer the complex parameter c to Mathematica under the name s. 

The parameters are stored in the form of rules in Mathematica, i.e. as name -> value. The 
first argument specifies the left-hand side of this rule. It need not be a symbol in the strict 
sense, but can be an arbitrary Mathematica expression. But note that in particular the 
underscore has a special meaning in Mathematica and may not be used in symbol names. 
The second argument is then the right-hand side of the rule and can be an arbitrary Fortran 
expression containing model parameters, kinematic variables, etc. 

The following example demonstrates the form of the parameter and data assignments. 
Shown are results of a differential cross-section for a 2 — > 2 reaction at one point in MSSM 
parameter space. Within the data the varied parameter is cos^, the scattering angle. 

Para[l] = { TB -> 1.5, MUE -> -1000., MSusy -> 1000., 
MAO -> 700. , M2 -> 100. } 
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Data[l] = { DataRow[{500. , -0.99}, 

{0 . 10592302458950732 , .016577997941111422} , 

{0., 0.}], 
DataRow [{500 . , -0 . 33} , 

{0 . 16495552191438356 , .014989931149150608} , 

{0., 0.}], 
DataRow [{500 . , . 33} , 

{0 . 2986891221231292 , . 015013326141014818} , 

{0., 0.}], 
DataRow [{500. , 0.99}, 

{0 . 5071238252157443 , . 012260927614082411} , 

{0., 0.}] } 

The DataRow [v , r , e] function has three arguments: 

• the independent kinematic variables (v = {\fs, cos 6*} above), 

• the cross-section or decay-rate results (r = {tree-level result, one-loop correction} 
above), and 

• the respective integration errors (e = {0, 0} above, as this example originates from 
the computation of a differential cross-section where no integration is performed). 

5 Using the Generated Mathematica Function 

To the Mathematica novice it may not be obvious how to use the function described above 
to analyse data, produce plots, etc. 

As an example, let us produce a contour plot of the cross-section in the M^o-tan/3 
plane. It is assumed that the function run has the two parameters MAO and TB in its 
argument list: 

Install ["run"] 

xs[sqrtS_, MA0_, TB_] := ( 
run[{sqrtS, sqrtS}, MAO, TB] ; 
Data[l] [[1,2]] 

) 

ContourPlot[xs[500, MAO, TB] , {MAO, 100, 500}, {TB, 5, 50}] 

The function xs runs the Fortran code and selects the data to plot. The first argument 
of run, {sqrtS, sqrtS}, instructs the Fortran code to compute the total cross-section for 
just one point in energy. We then select the first (and only) DataRow in the output and 
choose its second argument, the cross-section results: Data[l] [[1,2]]. 
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This example can be extended a little to produce a one- dimensional plot where e.g. for 
each value of tan (3 the minimum and maximum of the cross-section with respect to M^o 
is recorded: 



« Graphics 'FilledPlot' 

xsmin[sqrtS_, TB_] := 

FindMinimum[xs[sqrtS, MAO, TB] , {MAO, 100}] [[1]] 
xsmax [sqrtS_ , TB_] := 

-FindMinimum[-xs[sqrtS, MAO, TB] , {MAO, 100}] [[1]] 

FilledPlot [{xsmin [500, TB] , xsmax [500, TB]}, {TB, 5, 50}] 



6 Summary 

The new Mathematica interface makes it easy to turn Fortran code generated by FormCalc, 
thus far a stand-alone, command-line driven application, into a Mathematica function for 
numerically evaluating the cross-section. The generated function depends on any kinemat- 
ical or model parameter the user chooses. The Fortran code is accessed via the MathLink 
protocol. 

The advantages of such a function are obvious, as all of Mathematica' s sophisticated 
tools for visualization, numerical analysis, etc. can immediately be applied to it. Interfac- 
ing with Mathematica requires only minuscule changes to the parameter file run.F. The 
necessary changes in the build process are taken care of by the makefile. The interfacing 
code is available from FormCalc Version 5.1 onward. 

The author welcomes bug and performance reports, as well as suggestions for improve- 
ments at hahn@feynarts.de. 
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